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The collisionless interaction of an expanding high-energy plasma cloud with a magnetized back- 
ground plasma in the presence of a dipole magnetic field is examined in the framework of a 2D3V 
hybrid (kinetic ions and massless fiuid electrons) model. The retardation of the plasma cloud and 
' the dynamics of the perturbed electromagnetic fields and the background plasma are studied for 

, high Alfven-Mach numbers using the particle-in-cell method. It is shown that the plasma cloud 

' expands excluding the ambient magnetic field and the background plasma to form a diamagnetic 

cavity which is accompanied by the generation of a collisionless shock wave. The energy exchange 
between the plasma cloud and the background plasma is also studied and qualitative agreement 
with the analytical model suggested previously is obtained. 
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INTRODUCTION 



O I When a high-energy density plasma is expanded into a magnetized background plasma, the injected plasma expands 
and excludes the background plasma and magnetic field to form a diamagnetic cavity. This process has been observed 
in a number of experiments in space and the laboratory (see, e.g., Ref. [l| and references therein). The similar processes 
are responsible for the formation of the hehopause, where solar wind drags or convects magnetic fields from the Sun, 
^isj which causes the deflection of both the interstellar medium and cosmic rays. Other examples of naturally occurring 
f~i inflated or stretched magnetic fields include coronal mass ejections and the formation of the Earth's magnetotail 
Qh. 0, S Starting in the 1960s, many experiments of the collisionless phenomena of exploding plasmas in space and 
laboratory were performed or planned. Dense, high-energy plasmas can also be formed by laser-irradiation of a small 
target embedded in a gas in an external magnetic field ^5]- The background gas can either be pre-ionized or be 
ionized by the laser. Such experiments show that, in addition to expansion, the target plasma is also subject to fiute 
instabilities on its surface as it interacts with the background field and embedded plasma. The diamagnetic cavities 
'/"J • are expected to occur also in supernova explosions Q . The physics of the plasma expansion and evolution has been 

■ investigated in detail numerically, using a variety of techniques, such as full particle Q, hybrid (particle ions, fluid 
. electron) Q and magnetohydrodynamics (MHD) [Oj] codes. The properties of the expansion phase in the presence of 

■ either a stationary or a flowing background plasma [ll|, as well as the details of unstable modes on the surface 
of the expanding plasma [HI Hi US] ' have been well studied. 

Typical experimental parameters are such that the ions are collisionless but the situation for the electrons ranges 
from collisionless to collision-dominated. In this paper the dynamics of the expanding plasma cloud is studied for the 
completely collisionless regimes, which, in particular, are realized in many astrophysical processes (see, e.g., ^ and 
references therein). An example is the problem of the collisionless retardation of the supernovae remnants which was 
first formulated by Oort jT3| and subsequently analyzed in detail by Shklovskii p^ . 
' In recent years, the basic concept of the plasma expansion has been extended to consider a large magnetic bubble, 
which can be formed using a small magnetic coil and plasma source attached to a spacecraft, to efficiently inflate the 
bubble to a large cross-sectional area [1, [l^. A net force would be exerted on the spacecraft due to the deflection 
of the solar wind around the bubble. In this case, t he p lasma is continuously injected in the presence of a dipole 
or dipole-like magnetic field. Theoretical arguments [ij and calculations (see, e.g., Refs. [1, [H, [l^ and references 
therein) suggest that the magnetic field of the dipole can be expanded with the plasma, so that the magnitude of the 
field falls off much slower with distance r from the dipole, namely as , with s ~ 1, 2 in certain directions at least, 
rather than r^^ (for a bare dipole), allowing a large bubble to form. Laboratory experiments [lOiH^ have provided 
some evidence for the slow falloff of the field from the source. However, the nature of the plasma and magnetic field 



'Electronic address; |hrachya@irphe.am| 



2 



expansion in this configuration is not completely understood, and how it cornpares with the more common picture 
of diamagnetic cavity formation has not been addressed to date. References P, H, [l^ the papers cited therein 
conducted a systematic study in a two-dimensional (2D) geometry and only for initial phase of expansion 0, and 
in a three-dimensions (3D) hut mainly of the characteristics of the magnetic field inflation [l^. Special attention 
has been paid to the MHD analysis of the expanding plasma behavior in a dipolc field in a vacuum [2^ . i.e. in the 
absence of the background plasma. 

In this paper, in order to provide further insight into the underlying physics we have performed systematic particle- 
in-cell (PIC) 2D3V hybrid simulations of the plasma expansion in a magnetized background plasma in the presence 
of a dipole magnetic field. The basic parameters of the problem are introduced in Sec. HIl The numerical method and 
model used are discussed in Sec. IIIII We perform a number of simulations assuming hydrogen cloud and hydrogen 
background plasmas and the results are presented and discussed in Sec.|TVl Finally, we state the conclusions in Sec.lVl 

II. BASIC PARAMETERS FOR THE PLASMA EXPANSION 

The plasma expansion process is characterized by the magnetic (Rm) and hydrodynamic (Rg) retardation lengths. 
Rm is obtained by equating the initial kinetic energy Wq of an initi ally spherical plasma cloud to the energy of the 
magnetic field that it pushes out in expanding to the radius Rm [23|, i.e., i?„ = (6Wo/i?g )"'^/'^. Here Hq is the 
strength of the unperturbed magnetic field. As the plasma expands, it draws the background plasma into a combined 
motion. Along with this, the mass of expelled plasma increases. The radius of the sphere within which the mass of 
the plasma cloud and that of the background plasma drawn into the combined motion become equal is referred to 
as the hydrodynamic retardation length: Rg = (3M/47™hTOb)^/^, where nb and mt are the density and mass of the 
background plasma ions [lB| and M is the mass of the expanding plasma. The smaller of the lengths, Rm or Rg, 
determines the predominant mechanism for the retardation of the plasma cloud-magnetic or hydrodynamic. Of course 
Rm and Rg account for the ideal characteristics of the retardation process, therefore some high-energy particles may 

penetrate beyond the stopping length. The relation Rm/Rg ~ M^^^ , where Ma — u/va is the Alfven-Mach number 
[u - {W^/Mf/'^ is the initial expansion velocity of the plasma) and va = Hq / ^/Amri/rnh is the Alfven velocity in the 
background plasma, implies that for Ma ^ 1 the cloud loses energy as a result of the deformation and displacement 
of the magnetic field, while for Ma S> 1 the retardation is caused by the interaction with the background plasma 
[UjH^ (see also Ref. [l|). In this paper we consider only the second (hydrodynamic) regime of coUisionless expansion 
with Ma ^ 1. In this case the retardation length is given by ts — Rg- 

An analysis shows that the hydrodynamic retardation can only be ensured by a coUisionless laminar (or turbulent) 
mechanism [23| associated with the generation of vortical electric fields Ei in the front of the expanding plasma or by 
a coUisional mechanism owing to pairwise collisions of the cloud ions with the ions and electrons of the background 
plasma. The ions of the expanding plasma transfer their energy to the ions and electrons of the background plasma 
in multiple ion-ion or ion-electron Coulomb collisions with a mean free paths \u and Ag; ~ {me/mc)Xii, respectively 
[26l. [27j. Here TOe and TOc are the electron and cloud ion masses, respectively. In this paper we assume a completely 
coUisionless regimes when Xu 3> Xd > is- Strictly speaking, the hydrodynamic description is inapplicable under 
these conditions. However, as shown in Ref. [28| the hydrodynamic approach, when considering coUisionless plasma 
expansion into the magnetized, ionized medium, can give quite reasonable qualitative and even satisfactory quantita- 
tive results. The physical reason underlying this situation is the small ion cyclotron radius in comparison with a 
characteristic flow scale is, with the ion cyclotron radius playing the role of particle mean free path. 

The first group of the coUisionless interaction mechanisms are collective turbulent mechanisms (anomalous viscosity, 
anomalous resistivity) during the development of ion-ion or electron-ion beam instabilities [29l] . The condition for 
excitation of the ion-ion instability has the form u"^ v\ + 2c1 (see, e.g., Ref. [lOl) or < 1 + 2c1/v\, where 
Cs = \/T^Jrnb is the ion sound speed in the background plasma. In many typical situations va ^ and Ma ^ 2. 
Thus, according to the criteria mentioned above at high Alfven-Mach numbers Ma ^ 1 the retardation of the 
expanding plasma cannot be caused by the turbulent mechanisms. 

The second group is the coUisionless laminar retardation mechanism associated with the generation of vortical 
electric fields. It is known that the role of vortical electric fields becomes predominant as Ma increases, since 
Ei/Ep ~ 24], where Ep is the polarization electric field which arises as a result of the drop of the hydrodynamic 
and magnetic pressures at the front of the plasma cloud. A model for energy exchange between the cloud and 
the background plasma due to the combined effect of the gyromotion of the ions and the generation of vortical 
electric fields when Ma > 1 (so-called magnetolaminar mechanism (MLM)) has been proposed in Ref. Analytic 
solutions for the initial expansion phase, when only a vortical electric field Ei develops, showed that the fraction 
of energy transferred from the cloud to the background plasma is proportional to S = {Rg/aL)^ (MLM interaction 
parameter), where is the cyclotron radius of the cloud ions. 
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In the framework of the ideal MHD approximation for the description of the plasma expansion into a vacuum in 
the presence of dipole magnetic field, another energetic parameter k was defined which is given by k = Wq/Wm 
[sl ]. where Wm is the total magnetic energy of the dipole beyond the spherical radius [Wm = rd is the 

distance from the dipole to the center of the plasma cloud and p is the magnetic moment magnitude. Nikitin et 
al. introduced critical value Kc for a different plasma location [331 . the case when k < Kc, a substantial plasma 
retardation will occur in all directions from the cloud center location ("quasi-capture" mode), meanwhile the plasma 
will not be captured by an ambient magnetic field when k > Kc ("rupture" mode). The critical value varies between 
0.1 and 0.4. The latter value is realized when the plasma cloud is located at the dipole axis. In this paper we consider 
only the second regime of the plasma expansion. 



III. COLLISIONLESS HYBRID SIMULATION 



The hybrid model is used to study the coUisionlcss plasma expansion in a magnetic field. This model describes the 
ions by means of the velocity distribution function /(r, v, t), the electrons being considered hydrodynamically [3^.|3^. 
The typical scale on which the macroscopic parameters (the plasma density, magnetic field) vary is the ion cyclotron 
radius a^. The electron cyclotron radius a^e is much less than due to the small electron mass. Of course, the 
electron distribution function changes sharply on a small scale a^ej but we do not consider such details. If we are 
interested in the behavior of those macroscopic parameters that vary on a scale a^, it is sufficient to consider only 
the motion of the electron small-scale gyromotion centers. In this case one can consider the electrons as a fluid and 
describe their motions hydrodynamically. 

All of the estimates show that the plasma under the conditions in question is essentially quasineutral everywhere, 
although there are thin space charge regions. These regions are located near the plasma cloud boundary as well as 
near the coUisionless shock in the background plasma. The scales of these regions are essentially less than ol, and 
therefore, they are not considered here. Then the densities of ions and electrons rie in quasineutral plasma satisfy 
the condition — 2^71^, where Z^e is the charge of the ion from plasma species a. Let the cloud ions have the 
index a — c, whereas the index of the background ions is a = b. The ion distribution functions fa are governed by 
the Vlasov kinetic equation 

^+v.fi + ^fE+i[vxH]).f^^O. (1) 
ot or nia \ c J ov 

The electric and magnetic fields satisfy the Maxwell equations 

„ Annee , ^ „ ^ 1 9H _ , 

V X H V, - ve , V X E = — 2 

c c at 

in which the displacement current is omitted. We assume non-relativistic expansion velocity of the plasma cloud. 
Here is the ion mean velocity which is defined as 

Ue^'^ZaUa, Vi — ^ Z^nQ Vq , (3) 



(4) 



The mean electron velocity Ve satisfies the equation 

E = -i[veXH] (5) 
c 

which can be derived from the hydrodynamic equation of motion of electrons ignoring the inertial term (i.e. formally 
by taking the limit rrie — > 0). The electric field in Eq. (O is that required to keep the electrons and ions together 
ensuring the quasineutrality of the plasma. Thus, we do not resolve the expansion process on the Debye scale. One 
can interpret Eq. ^ as follows. The electromagnetic force acting on the electrons and ions does not depend on the 
particle mass; therefore light electrons are accelerated much more than ions. This would lead to charge separation and 
would violate the quasineutrality condition if the strong Coulomb interaction would not prevent this from occurring; 
i.e., an appropriate electric field arises to enforce quasineutrality. The mean force itself, which so strongly accelerates 
the electron gas, has to disappear or, more exactly, its electric component has to compensate the Lorentz force within 
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FIG. 1: Scheme of the simulation region showing location of the dipole and the initial position of the plasma cloud. 

the accuracy of the order of mg/mQ. This fact is expressed by the approximate relation ([5]). Equations 
constitute a closed system and are used to investigate the collisionless plasma dynamics in an external magnetic field. 

Note that in Eqs. (U)-® we have also neglected the terms associated with the finite conductivity (e.g. magnetic field 
diffusion and Joule heating). It can be shown that this approximation is easily satisfied in a wide range of parameters. 
In fact, for example, we can estimate the characteristic diffusion time for the magnetic field, Tm = Airaj^cr/c^ , where 
tJ is the conductivity of the plasma. Introducing the characteristic retardation time of the plasma cloud Tg = Rg/u 
we obtain Tm/Ts ~ M\{Xii/Rg) ^ 1. Here Xu is the mean free path for ion-ion Coulomb interaction. 

These approximations constitute the hybrid collisionless plasma model employed in this paper. It should be em- 
phasized that earlier simulations and experiments demonstrated the practical realizability of the magnetolaminar 
interaction model |2^] discussed briefly in Sec.|TTl In addition it can be shown that Eqs. H])-® of the hybrid model 
have universal nature. This universality can be demonstrated by normalizing Eqs. ([I])-® and all the quantities with 
the scales n;,o (densities), c/uopb (lengths), va (velocities) Hq (magnetic field), {va/c)Hq (electric field). In these units 
the time and the ionic distribution functions are scaled by the cyclotron period fi^ and nbo/v\, respectively. Here 
LOpb, fib and Ubo are the plasma and cyclotron frequencies and the unperturbed (initial) density of the background 
ions, respectively, va is the Alfven velocity in the background plasma. Then it is seen that Eqs. ([I])-® depend only 
on the dimensionless parameters Zc/Zb and Zcmb/Zbrric. 

Let us consider axially symmetric model for the dynamics of a point explosion that forms a cloud of dense plasma 
expanding into a magnetized low-density background plasma. For a simplicity we assume hydrogen cloud and back- 
ground plasmas with tUc = nib — ^Tip, Zc — Zb = 1, where rUp is the proton mass. At the initial time t = 0, an 
explosion occurs at the point p = z = in a cylindrical region ^ p ^ Lp, —L^ ^ z ^ Lz (see Fig. [1} filled with a 
homogeneous background plasma of density Ubo- The magnetic dipole with a moment p = Pz^z creates a magnetic 
field Ho(r) = Ho(p, z) and is placed outside the cylindrical region on the z-axis {zd < —Lz, = 0). Here p,z,ip 
are the cylindrical coordinates and p is the cylindrical radius vector. Similarly r, 9, ip are the spherical coordinates 
with the radius vector r, and 9 is the angle between r and the z-axis. We consider the case when the dipole moment 
is parallel to the z-axis and due to the symmetry we chose Pz > 0. In our simulations the dipole magnetic field is 
determined by the ratio rj = i?o(0, —Lz)/Ho{0,Lz) > 1 and the strength of the dipole magnetic field at the center of 
the simulation domain, i.e. Hq = Hqz{0) > 0. Using these two quantities it is easy to obtain the position Zd and the 
moment pz of the dipole 



which completely determine the dipole field for given length Lz ■ At fixed Hq the limits 77 — > 1 and 3> 1 correspond to 
nearly homogeneous and strongly inhomogeneous dipole magnetic fields in a cylindrical region, respectively. Through- 
out this paper the cyclotron frequency fib and the Alfven velocity va are determined in terms of the magnetic field 
Hq. Note that due to the axial symmetry all physical quantities are independent on the azimuthal angle ip (obviously 
this symmetry breaks down for arbitrary orientation of the magnetic moment p) . Here we neglect the (/^-derivatives 
in Eqs. ®-® and consider 2D3V model for our simulations. It should be also emphasized that the scheme outlined 
in Fig. [1] essentially differs from one adopted in Refs. 0, [l^ where the dipole is located inside the simulation domain 
and due to the strong interactions the ions may be trapped by the dipole. 

The explosion forms a cloud of dense plasma of radius rg containing Nc ions with a total kinetic energy Wq. At 
the initial time, the velocity of the cloud ions is distributed linearly along the radius, i.e. at r ^ vq 



Pz = ^Hq, Zd = -Lz 




(6) 



.1/3 _ 1 



Vr (r-,0) 



u(r) 



r 



(7) 




and Vr (r, 0) = at r > tq . Here is the maximal velocity of the ions in the cloud and is determined by the 
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TABLE I: The values of the parameters for the numerical simulation {Ma ~ 15, 5 = 1). The number of the ions (protons) in 
a cloud is iVc = 1.67 x lO^''. 





= 34.16 cm 


Rg 


= 34.16 cm 


Rm 


= 175.2 cm 


Ts = 


1.044 ^s 


Ho 


= 100 G 


Wq 


= 0.896 kJ 


Vm 


= 3.27 X 10^ cm/s 


V = 


100 


UbO 


= 10" cm-^ 


VA 


= 2.18 X 10^ 


cm/s Qb " 


= 9.58 X 10^ s'^ 


U!pb 


= 1.32 X 10^° s"^ 



initial energy Wo of the cloud, with i>„ (lOWo/STVcWp)^/^ = wo(5/3)i/2, where vq = {2Wo/Ncmpy^^ is the average 
velocity of the cloud ions. 

The equation of the energy conservation of the system is given by 

W = Wkin + 27r y pdp J dz i-nmvl + —j . (8) 

Here n = rie, and n is the total density of the ions in the cloud and background plasmas. The total energy of the 
system consists of the kinetic energies of the electron gas, the energy of the magnetic field, and the kinetic energy 
Wkin of the ions. Note that in Eq. ([8]) we have neglected the energy of the electric field as Ve -€1 c and E H. 

The system size employed in the computations is 2Lz = Q.Sa^^ and Lp = 4.90^ and uses 98 x 49 cells in z and p 
directions, respectively, where aL — Vm/^t is the ion cyclotron radius. Thus, the cell size in the p and z directions is 
h = hp = hz = O.laL- In order to correctly resolve the ion gyromotion the quantity h and the time step At are fraction 
of the ion cyclotron radius ol and the ion cyclotron period f^^^, respectively. In our simulations At — 0.01f2^^. The 
initial radius tq of the plasma cloud is considerably smaller than the step size h; that is, we can assume that at the 
initial time the cloud is concentrated at a point. At the initial time i = we take the background plasma at rest and 
assume that the background particles are distributed uniformly over the entire cylindrical region and the particles 
of the expanding plasma are distributed uniformly in the cloud. Explicitly the initial distribution functions of the 
background and cloud ions are given by /b(r,v,0) = ntoS{v) and /c(r,v, 0) = ncoS{v — u{r)^), respectively, where 
rico is the initial density of the cloud ions and u{r) is determined by Eq. ([7]). Note that the distribution function 
/c(r, V, 0) is non-zero only at ^ r ^ tq. Boundary conditions are one of the most important conditions to be satisfied 
in an electromagnetic problem, and they vary depending on the problem. In Ref. periodic boundary conditions 
are imposed on the fields and the particles. They only run the system a relatively short amount of time before any 
disturbances reach the edges of the system. In this paper, at the boundaries of the cylindrical region p = Lp and 
z = ±Lz, all quantities are specified by its unperturbed values, and on the z-axis (i.e. at p = 0) we assume the 
condition Ep^^p — Hp ^p — v^p^tp — which is a natural consequence of the symmetry of the model. The absorption 
boundary conditions are chosen for all the particles which means that any particle leaving the computation domain 
is assumed to have gotten lost. With these boundary conditions the calculations are continued until the time the 
perturbations reach the boundaries of the cylindrical region. 

The equations of motion for the ions are the equations of the characteristics of the Vlasov kinetic equation ([T|). 
Assuming hydrogen plasma we obtain 

E + i [v, x H] j . (9) 

These equations are solved by the particle- in-cell (PIC) method using Boris pusher algorithm [s^, [s^. Instead of 
the kinetic equation integration, the trajectories of a large number of the quasi-particles (each of these particles in 
turn consist of the large number of ions) are computed. In our simulations the total number of the quasi-particles 
is 97000 and 20168 quasi-particles are used for a cloud. For the background plasma we set 16 quasi-particles in 
the cell with close values of the cyclotron radii (because h <^ a^) which allows us to simulate ion kinetics in details. 
Further increasing the number of the quasi-particles does not essentially influence the simulation results. The Maxwell 
equations are solved using an explicit first order splitting scheme on staggered grid. A more detailed description of 
the mathematical model and the numerical simulation can be found, for example, in Refs. [s^ [35l[36|. 



IV. RESULTS 



In this section, results from several calculations are presented to describe the overall physics and to show the 
properties of the resulting structures. Numerical simulations have been performed for calculation of the energy, density 
and electromagnetic characteristics of the expansion of a hydrogen plasma into a uniform magnetized background 
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FIG. 2: (Color online) The time (in units of r^) variation of the kinetic energy (in units of Wa) of the plasma cloud and 
background plasma particles for Ma = 15.0, & = 1.0. a) The variation of the total energy of the background plasma (squares) 
and plasma cloud (circles). The dashed line is the theoretical predictions of Ref. [sJl (see the text for details), b) The time 
variation of the radial (Wp, circles), longitudinal {Wz, squares) and rotational (W^, triangles) energies of the plasma cloud, c) 
The same as in b) but for background plasma. 



hydrogen plasma at high Alfven-Mach numbers with Ma — Vm/vA — 15.0 and for a MLM interaction parameter 
6 = 1.0. The basic physical parameters of the simulations are shown in Table[Tl where = Rg/vm is the hydrodynamic 
retardation time. Note that Rg <C Rm and the retardation of the cloud is caused here by the interaction with the 
background plasma. The dipole is placed at Zd — —258.86 cm (see Fig. [1]). The initial radius of the plasma cloud 
is ro = O.Olai, which corresponds to the initial density rico = 10^" cm~^. The parameter k introduced in Sec. [Ill is 
K ~ 0.62 > Kc- Thus we expect that the ions will not be captured by a dipole magnetic field. Using the parameters 
introduced above, the simulation is run to t ~ St, when the perturbation reach to the boundaries of the cylindrical 
region and the boundary conditions adopted here become invalid. Due to the cylindrical symmetry the solutions of 
Eqs. dUl-ilS]) are obtained using the cylindrical coordinate system shown in Fig. [1] 

Figure [2] shows the time variation of the kinetic energies of the cloud and the background plasma. The energies are 
normalized to the initial energy of the plasma cloud Wq and the time is given in units of the retardation time . It 
is seen that the energy of the plasma cloud decreases with a time and at t ~ 5ts , about a half of the initial energy 
is transferred to the background plasma. Let us recall that the Coulomb collisions between particles are completely 
ignored here and thus as mentioned in Sec. |TT] the energy transfer is caused by the MLM interaction. In Ref. a 
simple analytical model was suggested for the energy transfer from the plasma cloud to the ambient plasma in the 
presence of a homogeneous magnetic field. Assuming the initial distribution functions /c(r,v, 0) = Nc6{r)fo{v) and 
/{,(r,v,0) = ri6(5(v) of the ions in the plasma cloud and the background plasma, respectively, the energy transfer in 
this model reads (see Ref. [sT] for details) 



AW{t) _ S 
Wo ~ 6 



where 



RgX 



2/3 



(10) 



(11) 



Here the velocity distribution function /o(w) is normalized according to V'(O) = 1 and the quantity K[t) is determined 
from the equation 



RgK{t) 



(12) 



Assuming that the strength of the homogeneous magnetic field is and the velocity distribution function /o(w) = 
(3/47rw^„)i?(w„i — v) (where H{z) is the Heaviside unit-step function) the prediction of this model is shown in Fig. 
(the dashed line). The transferred energy ((TU]) at i < and for a chosen /o(w) is then given by Aiy(t)/Wo — 
{S/6){t/Ts)^ and in the hmit t — > oo approaches AWoo/Wq = C6 (where C ~ 0.59) independently of the particles 
initial velocity distribution function in a cloud. It is seen that the theory agrees qualitatively with the simulation 
(squares) although it has been derived under assumption that AW{t) ^ Wq. The sum of the cloud and background 
plasma energies is the total kinetic energy Wkin(i) of the ions involved in the energy conservation relation, Eq. ([5]). 
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FIG. 3: (Color online) The perturbations of the electromagnetic fields and the densities of the plasma cloud and background 
plasma as functions of the spherical coordinate r (in cm) for t = 2rs, Ma ~ 15.0 and 5 = 1.0 and for some values of 9. The 
panels la-3a show the distributions of the p (solid lines), ip (dashed lines), and z (dotted lines) components of the electric field 
(in units of {va/c)Ho) at 6 — 7r/4, 9 — 7r/2 and 9 — 37r/4, respectively. The panels lb-3b show the same as la-3a but for 
the perturbation of the magnetic field AH (in units of Ho)- The components of the dipole magnetic field Hqz and Hop are 
also shown as thin solid lines. lc-3c show the densities (in units of nto) of the cloud (solid lines) and the background plasma 
(dashed lines) ai 6 — it/ 4, 9 = n/2 and 6 = 37r/4, respectively. 

This energy decreases with a time due to the energy transfer to the electrons (second term in Eq. ([5])) and the magnetic 
field (last term in Eq. ([H])). However the amount of the energy gained by the electrons and the magnetic field is small 
compared to Wkin- 

The time variation of the radial (Wp), longitudinal (Wz) and gyromotion (W^) components of the energies are 
shown in Fig. ^^,c. Using Eq. ([7]) it is easy to calculate the radial and longitudinal energies of the plasma cloud 
at t = which are given by Wp{0) = {2/3)Wq and M^z(O) = (l/3)Wo, respectively. At the initial time there is no 
gyromotion of the particles and W(jp(0) = 0. As shown in Fig.[2]3,c the expanding plasma cloud loses its radial energy 
which is transferred to the gyromotion energy of the cloud and the background plasma ions. In addition a large 
amount of the cloud kinetic energy is transferred to the radial energy of the background ions. It is also seen that the 
changes of the longitudinal energies of the ions are small. 

Figure [3] shows the electric (panels la-3a) and the perturbed magnetic (AH(r,t) = H(r,t) — Ho(r), panels lb~3b) 
fields strengths in units of {va/c)Ho and Hq, respectively, as well as the densities (panels lc-3c) of the cloud and the 
background plasma in units of ni,o as the functions of r (in cm) at t = 2ts, = 7r/4, = 7r/2 and 9 = 37r/4. Here 9 
is the angle between the spherical radius vector r and the z-axis. In the panels lb-3b of Fig. [3] the radial H^p and 
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FIG. 4: (Color online) The evolution of the boundary of the plasma cloud from t = O.IIts (central circle) to f ~ 5ts (external 
dashed line). The time interval between two successive lines is approximately At ~ OAts- 

longitudinal i/pz components of the dipole magnetic field are also shown (note that TJqv — 0) thin solid lines. As 
expected the density of the plasma cloud is strongly reduced due to the expansion. In addition the perturbations 
of the electromagnetic fields are larger at 6 ^ 7r/2 where the dipole magnetic field is stronger. At the initial stage 
of the expansion process {t <C Ts) the ions of the plasma cloud front propagates through a quiescent background 
plasma, and different plasmas are mixed. At this stage the background ions are accelerated mainly by the generated 
azimuthal electric field. Then, turning under the Lorentz force action, they acquire the radial velocity. At the later 
stage when the radial and longitudinal electric fields arise, the motionless background ions acquire a velocity due to 
all components of the field, with the azimuthal and radial fields being dominant. In general the motion of the ions in 
the electric and magnetic fields is very complicated. They will do a gyromotion along the magnetic field lines. But 
along with the gyromotion, the ions will also experience E x H drift and gradient drift. In contrast to the azimuthal 
electric field, the radial and longitudinal components of E are space-charge fields. This follows from Eq. ([5]) as the 
electromagnetic induction Eq. Q provides the azimuthal field generation rather than the radial and longitudinal ones. 
The electric field is generated in a region of the perturbation of the magnetic field, and everywhere it is perpendicular 
to the total field H. 

From Fig. [3] it is seen that the retardation of the plasma cloud is accompanied by the formation of a compressed 
plasma layer which moves together with the compressed electromagnetic field; i.e., a coUisionless shock wave is formed. 
The depth of the compressed layer is of the order of the ion cyclotron radius a^, in agreement with an estimate A ~ 
of the width A of the coUisionless shock front given in Ref. [s^l ■ The ions of the background plasma are pushed out 
of the expansion region; this leads to the formation of a plasma cavity. It correlates with the diamagnetic cavity, a 
region in which the total magnetic field H(r, t) is smaller than the unperturbed dipole field because it is squeezed 
out (Fig. [3]). Thus, there is essentially no electric field in the cavity. Due to the retardation of the ions of the plasma 
cloud front, a part of its mass forms a shell near the front boundary, see panels lc-3c of Fig.[3l However this shell is 
not isotropic and is more visible at 9 ^ 7r/2, i.e. in the region closer to the dipole. At the later time of the expansion 
t > 2ts the plasma shell becomes more and more pronounced with the formation of the cavity also in the plasma 
cloud with sharply (with the width ~ a^) distributed ions. 

The time evolution of the plasma cloud boundary is shown in Fig. [4l where the central circle and external dashed 
lines correspond to t = O.IIts and t ~ 5rs, respectively. As shown in Fig. 3] an initially symmetrical plasma cloud 
begins to be asymmetric in the pz plane at t > 2ts . On the other hand the shape of the cloud is symmetrical at all 
the stages in the plane perpendicular to the z-axis (not shown). The similar shape of the plasma cloud is observed 
in the experiments as well as in numerical simulations of the plasma expansion in a vacuum when the background 
plasma is absent [l|, [2^ . In this case the asymmetrical pattern may be caused only by the ambient magnetic field. 
At the initial stage, the plasma cloud expands isotropically as a free stream with high kinetic energy. Then, plasma 
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expanding downward begins to decelerate because the ambient dipole magnetic field is stronger at a lower region near 
the dipole (see, e.g., Figs. [T]and[l]). The kinetic beta Pc defined as the ratio of plasma cloud kinetic energy density 
to magnetic energy density reaches unity around this area first, and then strong interaction between plasma and the 
ambient field occurs. Thus, the diamagnetic motion of plasma is induced to generate the surface diamagnetic current 
asymmetrically, only on the lower plasma surface. In Rcf. 38] an analytical model was suggested to calculate the 
pressure of the dipole magnetic field on the surface of the plasma cloud. It was shown that at p = pzBz this pressure 
vanishes at = 0; tt and has its maximum at 9 = 0max (""/S < ^max < ti") which depends on the ratio R/\Z(i\. Here R 
is the radius of the plasma cloud. The maximal value 0max tends to tt with increasing R and shifts towards the value 
6'niax — 7r/2, when R ^ 0. Therefore the layer near 9 ~ ^max of the expanding plasma sphere will be mainly deformed 
by the external magnetic pressure. In the present context of the high-energy expansion /3c S> 1 even at the later 
stages of the expansion and the magnetic field cannot deform directly the cloud shape. In this case the deformation 
occurs due to the interaction with the background plasma with initially low kinetic f3b. Then the background plasma 
tends to change the shape to follow the dipole magnetic field lines. As shown in Fig. [2] the kinetic energy of the 
background ions and hence the parameter /3b increases with time which leads to the formation of the asymmetrical 
kinetic pressure of the background ions on the plasma cloud. 

V. CONCLUSIONS 

In this paper we have investigated the coUisionless expansion of the dense plasma cloud into magnetized background 
plasma in the presence of a dipole magnetic field. The 2D3V hybrid-PIC simulations with kinetic ions and massless 
fluid electrons have been employed. We have considered the case of high-energy super- Alfvcnic expansion with 
Ma ^ 1 when the expanding plasma is not captured by the ambient magnetic fleld. It is shown that in this 
parameter regime the retardation and the deformation of the plasma cloud is mainly caused by the interaction with 
background plasma which, however, can be strongly affected (at least at the initial stage) by a magnetic field. Our 
numerical results for the energy transfer from the plasma cloud to the background plasma have been compared with 
the theoretical predictions of Ref. [3l| and qualitatively good agreement has been found. For future applications 
the arbitrary orientation of the dipole (with respect to the z-axis, see Fig. [T]) should be considered. For instance, 
such a configuration has been realized in the experiment [2^ . In this case the physical quantities also depend on the 
azimuthal angle ip and the shape of the plasma cloud becomes asymmetric in the plane perpendicular to the z-axis. 

An additional item for further investigations is the development of a model which accounts the thermal effects 
for electrons which are completely neglected in the current study. These effects can be included by adding in the 
right hand side of Eq. ([5]) the term — (l/en)Vpe, where Pe is the thermal pressure of the electrons. Also along with 
Eqs. (P)-® an equation for the evolution of Pe should be considered (see, e.g., Ref. [1^). In this connection assuming 
the strong heating of the electrons in the shock front it is also desirable to investigate the physics allowing for charge 
separation effects and the resolution of the electron Debye length scale. We intend to address this and other issues in 
the context of the plasma expansion in a separate study. 
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